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Estimates made of column densities of H 2 at ■ OP indicate 
that H 2 is far more abundant than HI in the inner galaxy and is the 
key to an explanation of the y-ray observations. This is also 
reflected in the correlation of galactic longitude and latitude 
distributions of y-rays and molecular clouds. Particularly strong 
evidence is found from the galactic survey of CO emission at 
2.64 mm. Th« cosmic-ray distribution inferred from the calculations 
is not uniform but only weakly dependent on the total gas distribution 
in the inner galaxy. 

1. Introduction. Molecular hydrogen has long been suspected to be an important 
compo.ient of interstellar gas because it is the most stable low- temperature form of 
the most abundant element in the galaxy. H 2 is expected to be the predominant form 
of hydrogen in cool clouds of sufficient density. However, despite its obutidancu, it 
is difficult to measure its galactic distribution directly. Results from two pro- 
mising methods for indirectly studying the galactic distribution of II 2 are discussed 
here* viz.* recent galactic surveys of \00 MeV y-radiation and 2.6 mm radio line 
emission from the J • 1 O transition of CO molecules. To these sur 'eys, which 
reflect the extent and distribution of H 2 in the plane of galaxy* we will add 
corroborating information on the ^unount and latitude distribution of gas in the 
direction of the galactic center supplied by x-ray, optical* and infrared 
absorption measurrir.ents. 

2. The Recent SAS-2 Gamma-Ray Galactic Longitude Observations . Flchtel et al. 
(197S) have recently reported the results of a sky survey made of lOO MeV y-rays 
using a spark chamber aboard the SAS-2 satellite. The similar galactic longitude 
distribution observed by the OSO-3 detector was also noted to be distinctly uncor- 
related with the 21 cm distribution by Clark et al. (1970). The OSO-3 result 
implied an increase in cosr.’ic-rays* unseen gas or both in the inner galaxy* and it 
was suggested by Stacker (1969* 1971) and Stechei and Stecker (1970) that molecular 
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hydrogen unsnun in 21 cm surveyi could account for a large part of the y-ray 
enhancement in the inner galaxy. ^ model hae been proposed by Bignami and Fichtel 
(1974) and Bignami et al. (1975) based on producing large en)wincemonts at the 
locations of arms mapped by 21 cm surveys by posti:ulating higher gas densities than 
are seen in 21 cm and proportionately higher cosii.'ic-ray intensities in the arms. 
However, Burton et al. (1975) indicate that the arms outside of 25°fr(n tiio galactic 
center are optically thin in 21 cm and significantly higher amounts of HI in tJxs arms 
than t)K>se deduced in previous surveys appear to be ruled out. We, therefore, feel 
that the galactic y-ray observations can be better understood by using other 
observations in addition to 21 cut surveys to determine the role of H 2 clouds invisible 
in 21 cm omission. 

3. The Galactic CO Distribution and the Molecular Ring at *v>51tpc . A survey of the 
galactic longitude distribution of CO emission in the galactic plane has recently 
been made by Scoville and Solomon (1975). The importance of this survey in under- 
standing the distribution of H 2 in the galaxy lies in t)\e fact that CO is an 
excellent tracer of H 2 . 

Scoville and Solomon (1975) have used the velocity profile data obtained in 
their CO survey in conjunction with the Schmidt (1965) rotational model of t)ie 
galaxy to determine the mean distribution of CO in the galaxy as a function of 
galactocentric distance for distances greater than 2.0 )ipc. This distribution 
shows a broad pea)( with a maximum near 5 kpc which they have concluded indicates a 
ring of H 2 clouds in this region. The general form uf the imlerular cloud di.lriutL>n 
obtained by Scoville and Solomon has recently been confirmed in an independent CO 
survey by Burton et al. (1975). The connection between this feature and the y-ray 
omission ring at 'v5 )ipc (Puget and Stacker, 1974) led to the suggestion that the 
y-ray data also provide evidence for the molecular cloud ring near 5 kpc (Solomon 
and Stecker, 1974). Coincidently , there is also a similar distribution and peak 
in the giant Hll regions of the galaxy (Metzger 1970). This may be understood to 
>3e the effect of hot young stars in ng associations being formed out of dense 
molecular clouds in this ring. The formation of such a prominent molecular ring 
poses an intriguing problem for galactic structure theory. 

4. Molecular Hydrogen and Total Co.;.umn Densities in the Direction Various 

methods can be used to estimate the amourit of gas in the direction of the galactic 
center (Steckei. uc al. 1975). The results are summarized in Table 1. 


Table 1. Coluiiui Densities of Hydrogen at i * 0° Excluding the Galactic Nucleus 


( 

x 10"22) (cm”^) 

<Ng.C. ,0) 


rrom 21 cm 


Ko.& to 1.5 

Daltabuit and Meyer (1972) 

radio 

1 to 2 

Kerr and Westerhout (1965) 



<1.2 

Clark (1965) 

<2H„ > 
fr8m CO 


3 to 10 

Scoville and Solomon (1975) 

<2Nh2 + N„j> 


;i(ii.5±2) 

This work (Icr < I©) 

from SAS-2 

Y-ray flux 



<2Nh2 N,jj> 

from X-ray 

absorption 

6.5 -o 9 

OH 2 / 20 H 1 - 1*7 (Kaplan and Markin 
1973) as verified by the measure- 
ments of Crasemann et al. (1974). 

<2Nh 2 * N„j> 

from IR absorption 

5 to 7.5 

Ryter, et al. (1975) 
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Scovilla, Soloron and Jefferta (1974) antimate tha total maaa of tha molacular 
diak ttaar tha galactic cantar to ha 

4 X 10^ i ^ 

This ylalda an astimatad flux from decay In the range 

0.6x10"® iv ^^y,oc ^ 1.5x10"® cm"^ s"^ rad ^ 
which la only about one tenth of tha observed flux at l»0®, 

5. G a lactic Gamma Ray flux as a Function of Longitude . In parformlr j numerical 
calculations of tha longitude distribution of the golactic Y-ray flux from w® 
decay, we have used the survey of Scovllla and Solomon (1975) to obtain tha 
relative distribution of molecular hydrogen In the galaxy as a function of galacto- 
centric distance 'iV and have normalized to a total colunn density of 7x10^^ cm"^ 

In the direction of the galactic center using the X-ray and Infrared absorption 
measurements (Taole 1) which are consistent with, but presumably more accurate 
than the column densities deduced from the Y'tay or CO results. The contribution 
from atomic hydrogen was estimated based on the ntimbers given by Kerr and 
Westerhout (1965) and Westerhout (1970). 


For the purpose of the calculations to estimata tha effect of cosmic ray 
enhancements in the galaxy, It was assumed that such anhfncements may be correlated 
with the gas distribution so that 


J(u) _ kt iW ■»-2nH7(S) 


(7) 


22 o 

The results are shown in Figure 1 for “ 7 x 10 cm“* and a»0. 3 corresponding 

to an increase of about a factor of 2 in tne cosmic ray flux In the 5kpc region as 


indicated by studies of the supernova 
remnant distribution in the galaxy 
(Ilovaisky and Lequeux 1972, Kodaira 

1974) and the synchrotron radiation 
measurements (Daniel and Stephens 

1975) . This corresponds to a mean 
value for the total gas density of ~1 
atora/cm3 at 10 kpc of which MO per- 
cent is in molecular form. This 
agrees with the values given by Spltzer, 
et al. (1973) from measurements of 
..'otatlonal UV absorption lines of 1(2 
and those given by Jenkins and Savage 
(1974) for Lyman a lines of HI. In the 
5 kpc region, this corresponds to a 
volume-averaged total density of ''-5 
atoms/cm^ of which ''^O percent would be 
in molecular form. 

The results of the numerical cal- 
culations are shown in Figure 1 together 
with the flux distribution given by 
Fichtel et al. (1975) from the SAS-2 



Figure 1. Comparlso>' of the total 
calculated y-ray flux distribution 
with the SAS-2 longitude data. 
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ots«rvatlona for the half-plane from Cp to 180° over which the CO meaeuremente of 
Scovllle and Solomon (1975) can be applied. Effecta of tpcondary electron oroductiorv 
hremsetrahlung and Compton Interactions are Included. These results l'.dlcate that 
most of the observed Y-ray enhancement at low longitudes Is primarily a result of 
Increased gas density and that the molecular ring near 5 kpc plays an important 
role in accounting for this Increase. 

Cosmic rays are limited to a small variation over the galactic disk given 
by 0.?;^ n;J,0.5 as derived from the limits on the amount and distribution of gas 
implied by the observations discussed In section 4. We thus conclude that cosmic 
r ays cannot vary linearly with gas density over all segments of the galaxy nor 
do they appear to be uniformly distributed . 

The results Indicate that the character of the longitude distribution of 
galactic Y-radlatlon corresponds well with the overall density distribution In 
the galaxy Implied by the CO and 21 cm measurements and that this distribution has 
a broad maximum In the 5 to 6 kpc region. Higher frequency modulations by spiral 
arms do not appear to us to play a significant role In determining the galactic 
Y-ray diferlbutlon within the statistical errors and 5° resolution of the SAS-2 
data, at least for the half plane analyzed here using the CO data. In the oti.er 
half plane, 180°- 1 - 360°, three apparently sharp features have been identified 
by Plchtel et al. (1975) with the Scutum, Norma and "3 kpc" arm features designated 
in 21 cm surveys. The two Inner features at 330^-335° and 340°-345° may correspond 
to a bifurcation of the broad moleuclar ring near 5 kpc on the other half-plane 
of the galaxy; the feature at 310°-315° is perhaps somewhat more puzzling since 
this should be associated with a correspondingly strong feature at ^ “50° which 
appears to be absent in both the Y-ray and CO data There are, however, large 
error bars associated with the 310°-315° observation due to the fact that the 
SAS-2 spark-chamber telescope only viewed this direction obliquely. Also, Puget 
et al. (in p ap.) have Indicated that large corrections due t.o nearby features are 
necessary in the 310°-360° region. Future CO observations from the southern 
hemisphere could help Increase our understanding of the matter distribution In 
this region. 

Thus, the ~''actlc gas seems to have a large-scale superstructure modulated 
by spiral arm yvi turbatlons similar to that seen In M31 In 21 cm emission (Guibert 
1974, Emerson and in our own galaxy in nonthermal radio emi**®'- ■•'r:.— 197^ 

and It appears to be the superstructure which determines the character of the 
general central enhancement in the Y-ray longitude distribution. 

It would also appear that in both external galaxies and our own Galaxy, the 
density qradlent of the H 2 distribution is steeper than that of the HI distribution, 
resulting In a decrease in the ratio nH 2 /nHl with distance outward f“om the region 
of maximum density. 

The total contribution from bremsstrahlung and Compton interactions of both 
primary and secondary electrons to the galactic Y-ray flux in the central region 
Is of the order of 30 percent. In agreement with the estimates made by Stecker 
ut al. (1Q74) using the observed Y-ray energy spectrum obtained by SAS-2. 
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6. Gawioa-Ray Latitude Distribution and Aa^oclated Mna of Sight Rad<1«ntnq « In 
addition to the Y-ray longitude distribution measurements reported by Flchtel 
et al. (1975), a galactic latitude distribution was also obtained. This 
distribution shows an asymmetry with respect to the galactic plane with more 
flux coming from positive galactic latitudes in the case of moderate latitudes 
6*^ - lb I - 30^. This can be understood in the context of section 4 where it 
was pointed out that UV« X-ray infrared absorption measurements Indicate that the 
total column density of gas is proportional to the reddening in a given direction 
and that in the direction of highly reddened objects most of the hydrogen may well 
be in molecular fozm as expected in dense dust clouds and as is indicated by only 
a partial correlation with 21 cm emission. The data of Knapp and Kerr (197< ) in 
the region of the s)ty located between I ■ 345° and 3CP and between b ■ 10^ and 30P 
can be used to eva. uate the average reddening and HI column density at b ■ 20^. 

We obtain <Eg_y(20P)> ■ 0.26 mag and <Nj,(20°)* ■ 9.7 x 10^° H atoms cm“^. 

Postulating the same cosmic ray density as observed locally we find from the 
reddening and relation (2) a gamma-ray flux of 2.2 x 10”^ cm"^ g-l ^r”^ which com- 
pares favorably with the flux reported by Fichtel et al. (1975) observed at 
b ■ 20P of 2 to 2.5 X 10”^ x 10"^ cm~^ s“^ sr“^. ^wever, the flux estimated from 
the Ml-column density is only 1.4 x lO"^ cm”2 sec” sr"^ which is particularly 
significant in this case because at high galactic latitudes the 21 cm emission line 
is not optically thic)( (Knapp and Kerr 1974). Thus, the situation with regard to 
the latitude distribution of galactic gamma-rays is analogous to chat of the 
longitude distribution and can be txitter understood by talcing accou.nt of II 2 unseen 
in 2?. cm surveys of HI gas. 
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AS TRACERS OF lN^i*ERSTELLAR MATTER 

J. L. Pugot^' and C. Ryt«r^^ 

NASA/Goddard Spaca Flight Cantar 
Gr» tnbalt, Maryland U.8.A. 

G. Serra 

CESR, Toulousa« France 
and 

C. Blqnami^) 

Univarslta di Milano. Milano, Italy 

It ia ahown that tiia qaa to duat ratio in our galaxy 
ia tho same with a very good approximation even in rogiona 
of vory different dcnaitiee. The detailed comparison of 
the latitude distribution of the Y-ray intensity with 
reddening data shows that in the solar vicinity tits 
latter gives the best estimate of the total amount of 
interstellar gas. Thxs conclusion is used to compute 
the nearby contributions to tho Y-rays longitude distribution 
observed by SAS II and some implications for the large 
scale structure of tho galaxy are drawn. 

I. Introduction . It has become apparent in the last few years that the density 
of interstellar matter can deviate strongly from that of atomic hydrogen HI, 

as observed by its 21-cm transition line. An soon as the density reaches a 
few particles per cm^, hydrogen combines to form H 2 molecules, with the fraction 
of atomic hydrogen left being eventually as low as a few percent of molecular 
hydrogen in very dense clouds (liollcnbach et al. 1971; Heiles 1975). On a 
large scale, the effect is able to substantially alter tho picture of the 
gas distribution deduced solely from 21-cm observations. Fortunately, in dense 
regions where HI fails to be representative of the density, the column density 
of other elements may bo high enough to nuike several other absorption or emis- 
sion mochanibms easily observable. Among them are interstellar reddening, 
molecular emission or absorption, aiid Y-ray emission. Strictly speaking, each 
of these mechanisms is related to a well-defined physico-chemical species, 
i.e. interstellar reddening is to dust, centimeter or millimeter absorption 
or emission is due to molecules (other than H 2 ) and Y**ray emission is produced 
by the interaction of cosmic rays with nucleons, that is mostly with hydrogen, 
be It atomic, molecular, or ionized. 

II. Molecular Clouds . A typical molecular cloud can bo considered as formed 

of a dense core, with a column density of (5 lO) x 10^2 h atoms cm~2 and a 
size of one parsec > ^nd an extended halo, with a column density of (5 10) 

X 1q 2^ H atoms cm~2 and a diameter T>f 10 pc or so. Typical masses are 500 M^ 
for the core and 10^ M^ for the envelope. 

^NAS-NASA Research Associate 

2)nas/NA.SA Research Associate at the time of the conception of this work. 
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Th« tilBtrlbution of nearby (D < 1 pkc) nolacular clouds datactad to data 
la praaentod In figure 1. The survey might not be complete, but is probably 
fairly ropresentative. In figure 1, the nearby clouds located within 'V'loP frosi 
the galactic plane are represented by small bars as « function of longitude. In 
tlte lower part of the figure, the profiles of interstellar reddening in a 
b • ‘vicy° strip for distances 0.5 kpc (solid line) and 1 kpc (dotted line) of 
the sun, are reproaentod as they can be deduced from the compilation of Pits- 
gerald (19t,8) . The correlation between the color excess and the presence of 
clouds is ax^parent, and it can be concluded that molecular clouds are detected 
wlxjrevcr the color excess reaches Eb.y % 0.4 mag. According to eq. 1(3), the 
column density is thus Nh " 3 x lO^I H atom cm**. Interstellar reddening thus 
appears to be a faithful 
indiuutur of the preseni'c 
of molecular clouds, which, 
in turn accommodate roost 
of the interstellar roatter 
evo 1 in moderately • 

Ob' .cured regions 
(Hollonbach et al., 1971). 

Figure 1. Profxle of the average interstellar 
reddening for distances of 0.5 kpc (solid 
line) and 1 kpc (dashed line) of the sun. 



III. The Gas to Dust Ratio . Interstellar matter is made of hydrogen, mixed 
witit a few percent (by mass) of elements with A - 12, a fraction of them being 

condensed in dust. The dominant gaseous component is obviously hydrogen, but 
a substantial fraction of it can be either ionised, cr molecular, and thus es- 
capes 21-cro detection. The only attempt to include any form of hydrogen, a 
fact which may account for the large scatter in the gas-to-dust ratios found in 
the literature, has been made by Jenkins and Savage (1974). In their ..urvey, 
they obtained a value N^jj/Eq.v ■ 5.5 x 10^^ H atoms cm“2 mag”^, in agreement 
with most of the other measurements. They estimated that molecular and ionized 
hydrogen contribute 2 x 10^^ H atoms cm~^ mag*^ to reach their final value, 
which applic; to a region within less than '•'1 kpc of the sun. 

Ga.s-to-dust ratios can also be deduced in molecular clouds, where most of 
the hydrogen is molecular, using the theory of excitation of molecules. 

Ktltncr (1973) notes that in a cloud located in the Taurus complex, the formal- 
dehyde molecule K 2 CO observed in absorption correlates well with blue and red 
extinction. From the number ht quotes, one can deduce Nh/Eq.v ^ (6. 7) x 10^^ 
cm"2 where refers to the number of H atoms, that is Njx ■ Ziu,. More 

accur.itely, Encrenaz et al. (1975) find for carbon monoxide (CO) in emission 
the ratios Nco/^ " 1.6 x 10^^ molecules cm"2 mag"I and 2N||2/Nqq ■ 1.3 x lO^, 
as scon as Ay > 2. With Ay ■ 3.2 Eo-y, we obtain Nj^/E 3 _y ■ 6.6 x lO^i H atom 
cro~2 mag”l. Finally, the x-ray absorption, which mostly reflect the column 
density of M-eloments (C, N, O, Ne, ...) in the line of sight, was compared to 
the interstellar reddening of supernova remnants (Ryter et al. 1975). Postu- 
lating universal abundances as is usual when dealing with interstellar x-ray 
absorption, the value Nji/Es-y ■ (6.8 - 1.2) x 10^^ H atom cm~^ mag. has been 
found. We adopt the value 

Nh/Eq-h ■ 7 X lO^l H atoms cm”^ mag“I (1) 

i 
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IV. Cawm-My Ififion . T2i« Y~ray iMckgrouiu! !■ genor«lly b«liov«d 

to b« du« it ch« d«cay of neutral pions fomed in tha interaction of cosmic rays 
with Intarscallar matter. The distribution of the intensity thus reflects the 
product of tha density of interstellar matter and that of cosbic raya. A sodel 
uf tho distribution of matter and cossiio rays, based on the dyn.usics of tha 
Galaxy, iias Uien proposed by Bignami at al. (1975) to account for tho observed 
longitude distribution (/ichtal at al. 1975) i following a more empirical ap- 
proach, Stocker at al. (1975) suggested that the gamma-ray profile closely 
Bhitcltos the molecular hydrogen distribution deduced from 00 n,' lervations 
(Scovilla and Soloson 1975), at leant at positive longitude v' ire a molecular 
emission survey exists. Further refinements have been brough by Stacker (1975), 
who showed t)iat the required large scale cosmic denrity varia ion is best ro- 
pxaserited by taking it proportional to the supernova distribu ion in the Galaxy. 


All these treatments suffer from neglecting tho effect * I the olose-by 
interstellar medium, which is very patchy, and appears to ir roduce features 
which should be easily detectable with present Y-ray sansit vity. The aim of 
this section is to talce into account the local production r .' Y~rays in inter- 
stellar upace, that is within one kpc or so, and to avalua .e the consequences 
on the large scale distributio... 


The >-ray intensity ly associated with a line of s'ght is given by 

^ 10"26 Ny (ph cn-2 ,-l f-l) (2) 

using the production rate per hydrogen atom for a o'/smic-ray density equal to 
that observed in the solar vicinity (Stacker 1973). Ufinq Eq. 1, this relation 
can be expressed as 

ly - 7.25 X lo"® (ph cm"^ a"^ Br”^- mag”') (3) 


In order to case comparisons with otho; data, the ganwia-ray intensity can be 
expressed through inverting eq. 2 and writing Ny instead of N|{, i.e. 


Ny - lo26 atoms 

Figure 2. a upper part. latitude 
distribution of y-ray intensity 
(solid line), total hydrogen column 
density deduced from equation 1 
(dashed line), and ill c'Olumn density 
(crosses) . Data are averaged over 
fyCP in the direction of the 
galactic anticenter. *i 

b. lower part. The corrcsiK>nding 
distributions for tho inner ‘ 

Galaxy. Data are averaged over 
6CP in the direction of the 
galactic center. 
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TIm lAtitudo variation of th« iotanaity Xy AV«r« 9 «d ovar 60° in lonqitu<la, haa 
boon givan by richtal at al. (197S) for tha cantar and anti-cantar of tha 
Galaxy. Tha data ara raproducad in figura 2, a and b, along with tha total hy- 
drogen column dvnsity« daducad from an intarstallar raddoning compilation 
(ritzgorald 1968) through aquation 1, and avaragod in tha aama longitudr t vigai 
thu corresponding column dansitioa obtainad from 21-c« omi8sion« Nk|« tAV’>.'*>. from 
f'lltabult and Meyar (1970) ara also raproaonted. Inapaction of t)>a figura sug- 
goata tlio following commantai at latitudes larger than 15° w)icra tha ragions 
ubser/ad ara oxpactod to contain cosmic rays with tha same intensity as t)ie 
intensity obsarved in t)M> solar vicinity^ thara is a vary good agroomant bettraen 
the column density derived from rad laning and Y~ray data. Tha atomic hydrogen 
seen in 21 cm is seen to be only a fraction of t)M total density wliete molecular 
clouds are )cnown to be important! et low latitudes in the direction of tlw inner 
galaxy and in regions connected with Gould's Iwlt. The difference between t)te 
column density derived from reddening data at low latitude* show that the cosmic 
ray intensity must increase in the inner galaxy by about a factor of 2 on the 
average (a similar conclusion )ias been reached by Stocker at al. 1975 and 
Dlgnami et al. 1975) and decrease outside the solar circle as already noted by 
Dodds 'et al . 1975. Finally, it must be noted that the n° production in the 
Galactic plane at high latitude* contributes a flux Zy ■ 4 x 10"® ph cm"^ s"^sr"^ 
at t)u* (Xkles, which cannot be neglected when assersing the intencity of an 
extragalactic isotropic background. 

Thu cot elusion which is new and relevant to our purpose here is that the 
prop»'' 'tttality lM.‘twoen the Y~ray intensity ly and the interstellar reddening 
•xpresnod by equation 3, is %m11 confirmed by the presently available 
1 ’ ..ucio distribution of Y~ray wiicro the cosmic ray density is l>elieved to be 
uniiurm. We are thus confident t)iat it can be used to take into account tha 
srarll scale distribution of the closeby interstellar medium and deduce its 
contribution to the Y-ray distribution. 

Thu patchiness of closeby features along with the finite angula\ resolution 
of Y~ray detuctors may to soinu extent mimic t)ie presence of sources, which Yiavm 
first boon associated with molecular clouds by Black and Fazio (1973). It can 
also bo noteo that the column density of the core of dense molecular clouds is 
compar.ible t'' ''he column density to tho Galactic Center, and they would stand up 
as brigltt (ly % lo"^ ph cm"^ s""^ sr"i) if tho angular resolution of the 

detector wore high enough, tliat i-} better than 1/2°. The likeli)>ood of the 
existence uf V'^ay source-like excesses has been recently discussed by Serra and 
Niul (1975). 

1‘iuuru 3. Y"f*y distribution as 
■I function of galactic longitude, 

KiiinuiMHl up in ,i b ■ ± 10° strip, 
fhitcr corilour: original SAS II 

data from Kichtol et al. (1975). 
twack area: pre<Ucted contri- 

bution from the Intcrotcllar 
nwdiuiii within 0.5 kpc of t)x} sun, 
a*i drvjucod from interstellar 
reddening and equation 1. 

Sitadud area: contribution of 

the interstellar medium bet%#ecn 
0.5 and 1 kpc of the sun. 
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The roddeni'iq data in the atrip Ab t 10°« ahown in figure lb, have been 
used to predict the contribution c£ nearby features tc the Y-ray intensity ly* 
through equation 3. Tlia contribution of matter uithin O.S kpc of the sun (black 
area) ar.a from 0.5 to 1 kpc (dashed area) are represented in Figure 3 superimposed 
on the original profile obtained by Fichtel et al. (1975) from the ShS II data. 
Simiiaiily, in Figure 4, the longitude distribution computed by Stacker 
(1975) from CO data (solid line) and the sa>4e after correction for local featurus 
(dashed line) are mproduccd. The fit is obviously improved, and an analysis of 
the details of the Y-ray structure is quite appe>ling. Note that the corrections 
generally go in the right direction, in spite of a trend to bo slightly too weak. 

Figure 4. The ganma-ray distribution 
at positive galactic longitude as 
obtained by Stocker (1975) 

(soli! line), as corrected for nearly 
contributions (this work, dashed line) 
and the carbon monoxide distribution 
(Scoville and Solomon 1975). Nolo the 
improvement of the agreement between 
CO data and Y-ray distribution when 
corrected for nearby contributions. 


The blank area which is left in Fig. 3 now represents t)ie Y-ray emission 
profile produced by the largo scale features, allowing a discussion of galactic 
structure. A peak at the Galactic Center, with a width of order 15° (or ''-2 kpc), 
clearly stands up and may be associated with inverse Compton emission (Dodds 
et al. 1975). At f’^30°, a well-defined maximum becomes evident and exactly 
coii.b-sponds to the peak in the longitude distribution of CO (Scoville and Solomon 
1975, Burton et al. 1975). It a. so fits to the profile computed on the basis of 
a detailed model of the Galaxy (Bignami et al. 1975), and corresponds to the 
double feature around *• *'• 335". However, at f 30 ° both arms are too cloue to be 
resolved with a 5° resolution. It is interesting to note that the distribution 
of giant .III regions (Mezger 1970, Wilson et al. 1970) exhibit a very simila- 
overall property but with a more detailed structure. Finally, in the region 
no°<f <140°, the contribution from the outer Galaxy (beyona 1 *'< 2 kpc of the sun) 
ooeins to be very weak, a conclusion consistent with a sharp roll-off of the cos- 
mic ray density at R > 12 kpc. The existence of sharp features in the Y-ray 
profile indicates that matter have to be located in narrow lanes at the inner 
edge of the spiral arms defined from 21 cm measurements. Similarly, narrow dust 
lanes (obviously associated to molecular clouds in our Galaxv) are observed along 
arms in external galaxies (Lynds 1975) . The radial distribution of matter de- 
duced through unfolding of the longitude distributions implies a rapid roll off 
of the density of the lanes along the spiral arms as they wind out starting at 
M kpc. Such a picture is in agreement with the model developed by Ro)serts 
(1974) for spiral structure. 
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ON THE ORIrtlN OF COSMIC RAYS 


F. W. St«ckor 

TTmoretical Studies Croup, NASA/Goddard Space Plight Center 
Groenhelt, Maryland 20771 

Using recent surveys of molecular clouds arid y-rays in the 
galaxy, it has bacone possible to determine the distribution 
of 1 to 10 GeV cosmic-ray nucleons in the galaxy. This 
distribution appears to be identical to “he supernova 
r>mnant distribution to within experimental error, providing 
strong stipport for the hypothesis that supernovae produce 
most of the observed cosmic rays. This distribution resembles 
►hat of CD associations of average ago ''<30 million years 
suggesting that cosmic rays are produced by iK>pulation I 
objects about 30 million years after their birth. 

The problem of the origin of cosmic rays has been the central problem of 
high energy astrophysics for over a generation (see, e.g., Shapiro 196?, Glnrhurg 
and Syrovatskil 1904, Rosen 1969, and Precher and Fnirbidge 1972). 

Shortly after the discovery of the 3K microwave blackbody background radiation, 
it was first noted by Fazio, Stecker and Wright (1966) that such i.adiation precluded 
the existence of cosmic-ray electrons c ’tslde the galaxy with the same intensity 
as that observed locally. They noted tnat such electrons would produce more 100 
MeV Y~rays than observed if their distribution extended more than '^/30 kpc from 
the earth. Thus, it became apparent that cosmic-ray electrons were of galactic 
origin. 

It was further noted by Greisen (1966) and 7.atsepin and Kuz'min (1966) chat 
the blackbody background radiation would Interact with ultrahigh energy cosmic rays 
of extragalactic origin to produce a cutoff in their energy spectrum. The lack of 
an observed cutoff allowed one to rule out the universal origin hypothesis for 
ultrahigh energy cosmic rays and to place limits on the extent of their source 
region as being within 300 Mpe (Stecker 1968). Galactic origin cf even ultrahigh 
energy cosmic rays has been advocated (Stecker 1971, Syrovatskil 1971) and support 
of this hypothosi-i has been recently provided by indications of the anisotropy 
of ultrahigh energy cosmic rays (Krasilnokov et al. 1974, Hillas and Ouldrldge 1975). 

Attempts to place limits on the extragalactic cosmic-ray flux using the y-rny 
bacivgrourd observations have not been conclusive. They still allow a universal 
cosmic-ray nucleon flux provided the mean intergalactic gas density njQ < 10"^cm"^ 

and allow origin within the local supercli ster even if n^^ “ 10“^cm“^ (Stecker 1975) . 

Thus, the discussion of the extragalactic versus galactic origin hypothesis has 
continued down “o the present (Brecher and Burbidge 1972, Ginzburg 1974). 

It has long bean realized that observations of galactic y-xays .uuld provide 
Important information for resolving this problem, but until now it has not been 
possible because of insufficient y-ray data and an incomplete knowledge of the 
amount and distribution of an important component of the interstellar gas, viz., 
molecular hydrogen. Recent observations of the large-scale galactic distribution 
of Y“*^*diation (Fichtei et al. 1975) and molecular clouds (Scoville and Solomon 1975) 
have now made it possible to investigate the large-scale distribution of galactic 
cosmic rays. Using the new observations, Stecker et al. (1975) have determined 


that th« coamlc-ray distribution in the galaxy is not uniform as would be indicated 
by the extragalactic origin hypothesis. These results indicated that there is a 
weak correlation of the cosmic-ray flux with gas density* (mostly II 2 clouds) in the 
inner part of the galaxy. Also apparent was a falloff of the cosmic-ray flux in 
the outer galaxy (also found by Dodds, et al. 1975). The striking similarity 
between the cosmic-ray distribution deduced by Stecker et al. (1975) and the 
supernova distribution in the galaxy (Ilovaisky and Lequoux 1972, Kodaira 1975) 
provides new evidence that supernovae produce the bulk of t)ie cosmic-ray flux. 


The galactic y-raya are pri- 
marily the result of the decay of 
n®-mesons produced in cosmic-ray 
interactions with interstellar gas 
(Knlffon et al. 1970). Their flux 
is tlierofore proportional to the 
product of gas density and cosmic- 
ray Intensity Integrated along the 
line-of-sight within the solid 
angle subtended by the y-r&y 
telescope. If a cosmic-ray flux 
distribution is assumed and Compton 
and bremsstrahlung y-raya are also 
included in the calculation (a 30 
percent correction at most), one 
c^m calculate the flux expected to 
be observed by the SAS-2 yi&y 
telescope of Flchtel et al. (1975) 
integrated over +10° in galactic 
latitude and averaged over 5° 
longitude. This can only be done 
over the half of the galaxy for 
which the molecular cloud distribu- 
tion has been determined (Scoville 
and Solomon 1975) . The details of 
this calculation are given by 
Stecker u*- -.1. (1?75) . Tht y found 
that a unitorm cosmic-ray flux 
distributicu leads to a y-ray flux 
which is a factor of *'<2 too high 
compared to the observations in the 
anticenter direction and which is 
too low in the direction of the 
galactic center. 

Using the seune methods, a cal- 
culation of the y-ray flux distribu- 
tion can also be made under the 
as3um{Jtion that the cosmic-ray 
distribution is proportional to the 
supernova remnant distribution in the 
galaxy, as would be expected if (1) 
supernovae are the principal source 
of cosmic r-ys, and (2) cosmic rays 



P'gure 1, (a) (top) Galactic 
distributions of cosmic-ray intensity 
using the supernova remnant distribution 
of Kodaira (1974) and total gas density 
deduced by Stecker et al. (1975) . 

(b) (bottom) Distribution 
of ir® production rate In the galaxy 
based on (a) . 
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diffuse only a few hundred parsecs 
before leaking out of the galactic 
disk (Joklpll and Parker 1969, 

Pamaty, Poames and Lingunfelter 
1970) . 

I'eing the Buj.jerr'ova distribu- 
tion obtained by K^odalra (1974) an 
roprenentatlve of the oalactlc 
cosmic-ray flux distribution 

(normalized so ■ 1 at 10 kpc), 

wH 0 

the longitude distribution of 
Y-rays as would be observed by the 
SAS-2 telescope has been calculated. 

Figure la shows the cosmic-r ly and 
total matter distributions used 
(see Stecker ot al) . Figure lb 
shows the calculated Y"ray »!tr.ls- 
slvity from ir°-docay. The results 
are shown by the histogram in 
Figure ? along with the data 
actually obtained by Flchtel et al. 

(1975), Indicated by tl>e vertical 
lines. The calculated distribution 
is in remarkable agreement with the 
data, providing strong support for 
t)>e hypothesis that supernovao 
P» oduce most of the observed 
cosmic rays. 

The supernova remnant distribution used here can be correlated with the 
distributions of various other population I objects to estimate the average age of 
the remnants. For this purpose, we turn to the detailed discussion of the 
correlation of gas and OB associations made recently regarding M31 (Emerson 1974) . 
These results show that the surface densities of HII regions and atomic hydrogen 

2 23 

(HI) are related by “(o^^,) ’ . For our own galaxy, due in part to optical 

depth effects, the HI distribution obtained from 21 cm observations may be to some 
extent misleading in the inner galaxy (Scoville and Solomon 1975, Stocker et al . 
1975). However, using the data of Merger (1970) we find a correlation between the 
supernova remnant distribution and the HII region distribution of the form 

0 “(o Assuming, as for M31, o “(o )^*^^, it then foliows that 

&N rXX HXX HX 

c_„ *(o„_)®*®^. If the correlation between OB associations and HI gas is also 

expressed in the form c «(Ojj^)**', it ic found that the older the association, the 
smaller the value of m (Emerson 1974) . The decrease of this correlation with age 
is generally attributed to a spreading out of t)>e stars in the OB associations with 
time, all the stars in the association having been spawned from the seune cloud 
complex. The correlation defined by m ■ 0.89 is close to the mean correlation 
with gas of all OB associations and consistent with an age of 3 x 10^ yr. Thus, 
it appears that the sources of galactic cosmic rays are young galactic {copulation 

1 objects, most likely to be supvernovae. 



Figure 2. Calculated longitude 
distribution of galactic y-rays 
under the supernova origin hypothesis 
Oiistogram) compared with the 
observations of Fichtel et al. 

(1975) (verticjil lines) . 
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ThR IWTRRACrrONS OP yLTRAHIGH RNRIWY COSMIC RAY WCUSl VITM T!rrEP.OAI.ACTIC 
PHO'roN PIETJKi REUVTION BKTWEEN M?:AN ATOMIC MASS, AGE, AJfO ORIGIN 

J. L. Pugat* and P. N. Sneekar 

*!^ot«tio«I Studies Group, NASA/Godd&rd Space Plight Center 
Craenbelt, Maryland 20771 

v#e present noma results of detailed Monte Carlo calculations of 
the interactions bet%raen ultrahigh energy cosmxc-ray nuclei and 
intergalactlc radiation fields using improved estimations of 
these radiation fields and recent experimental determinations 
of cross sectlona for single and multiple photodisintegrations as 
a function of energy for 1 < A < 56. Starting from Fe, disintegra- 
tion product distributions are determined as a function of time. 

Som'4 of these results and their possible astrophyslcel implications 
are discussed here. 

1. Introduction. Of the various photon-nucleus Interactions possible between 
Icw-enargy intergalactlc photon fields and ultrahigh energy cosmic-ray nuclei, 
the only two of astrophysical significance for attenuation by energy loss In the 
lOl0-lo5l eV range are pair production and photodisintegration. The physics of 
pair production is well determined, and an excellent discussion of the relevrj^t 
cross section formulas has been given by Blumenthal (1970) which we have 
adopted in our calculations. The photodls integration process is not as easy 
to treat and most previous cosmte-ray papers have used only simple estimations 
of the crois section for thin process (CeraBlmovu and Pozontal * , 10C2» Grnisen, 
1966i Zatsepin and Kuz'min, 1966> Berezins)cii and Zatsepin, 1971| and Tkaezyk 
ct al. , In press) . 

Stacker (1969) used a detailed study (Gorbunov 1968) of the breakup of 
^Ha to determine the lifetime of ^Ue and ^®Fe against one nucleon loss by 
interactions with tha universal blackbody radiation. Tlie present paper presents 
some of the results of a new and much more detailed calculation of Intergaluctic 
photodisintegration of ultrahigh energy nuclei with the following Inprovemants i 

1) Use of detailed en^lrlcally determined cross section data as a function 
of energy for all nuclei with 1 < A £ 56. (Heavier nuclei are considered too 
rare to be of significance in the present context). 

2) Use of cross sectl«>n data for 1 £ A £ 56 for multlnueleon diaintegrat ons 
as well as single nucleon losses for all ems energies up to ‘'*140 MeV (the w 
production threshold) . 

3) Use of new estimates of the intergalactlc Infrared radiation field 
whose effect is Included along with l.iat of tha blackbody radiation. 

4) Inclusion of energy loss by pair production as well as photodisintegra- 
tlon in the calculation. 


•NAS/NRC Postdoctoral Resident Research Associate. Permanent Address i 
Observatolre de Meudon, Meudon,, France. 
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.>) U*« of numerical Mo"ta Carlo techniques in order to take proper account 

of multinucleon loss channels to determine the probability distribution of photo- 
disintegration histories of nuclei as a function of time. 


2. Calculations. The rate for photodisintegration of a nucleus of mass A with the 
release of N nucleons is given by (Stacker 1969) 


A,N 


- 1 

2 A 


r 

o 


dc 


nlc) 

c* 


/^"^'de'fc. (€•) 

0 


( 1 ) 


where n(e) is the photon density of the background radiation with energy e in the 
observer's system. 


While one and two nucleon losses are dominated, for the most part, by the giant 
resonance in the photonuclaon cross section up to a cms energy c' of -30 MeV. cross 
sections for N > 2 become important for A > 10 and c' > 30 MeV. For A -9. the 
disintegration is dominated by the reaction 


y * ^Be -► 2a n 


( 2 ) 


The background radiation is assumed to consist of three components! 


1) the blackbody radiation at T • 2.7 K» 

2) a dilute optical blackbody field T- 5000 K with a dilution factor 
of 1.2 lO'^Sf and 


3) a power-law infrared radiation field estimated from the contributions 
of extragalactic strong infrared sources such as quasars and Seyfert galaxies. 

Tv'o power-law forma were considered as probable upper and lower limits for the 
extragalactic infrared radiation based on cosmological models involving the 
evolution of extragalactic infrared sources. These power-law photon spectra. 
wh..ch were to extend from the optical graybody field to the microwave 

blackbody field, are given byi 

High Infrared! n(c) ■ 1.1 x 10“^ c 
(HIR) ® 

-3 -2 

Low Infrared! n(e) ■ 2.6 x 10 c „ 

(LIR) 

both taken over the range 2 x 10”^ S c S 0.83 eV. 

3. Results. Figure 1 shows the calculated loss time (inverse loss rate) for 
^*^Fe. '^Be and ^He as a function of Lorentz factor y for the HIR (solid line) and 
LIR (broken line) cases. The effects of the differing shapes of o (c ' ) can be 
seen as well as the significant difference in rates between the HIR and LIR cases 
(almost a factor of 10 for 56pa with y 3 x 10^) . Also shown are some partial 
channels for these nuclei so that the effect of multinucleon channels on the total 
Fe lifetime can be seen as well as the shape of the 2 nucleon contribution to the 
lifetime of He. 


(cm"^ eV ^) 


I 


I 
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Ttvi calculatft^ hlatory of a bean 
of nuolol atartin? at two 9 lvan 
•nartrlaa r.^ • AM,, ^ for tlia HIR 
caaa are stiown In Figure 2. Shown 
for compariaon are the lifetimes for 
Fe againat one nucleon loaa at each 
Y^« Xt la ohvloua that the llfetlmea 
for complete disintegration are 
conalderahly higher (a factor of 
^200 for the . .R caae) . 

The gap in the dlatribution for 

5 < A < 8 la caused by reaction (2) . 

The Monte Carlo calculations yielded 

the mean mans A and spread in mass 

(defined an a standard deviation 

from the mean) of the products of Fe 

disintegration as a function of time. 

Some results are shown in figure 3 

for three values of y . 

o 

The shaded areas for two of 
the curves illustrate the standard 
deviation spread in mass values. 

■^ese three curves illustrate the 
time evolution of Fa nuclei for the 
three basic energy regions where the 
physics is somewhat different. In 

the energy range y 2 5 x lO^, the nucleon loss timea for 

energy loss is dom?nated by photo- ^ including the effect 

disintegration reaction with the 2.7 nultinucleon channels. TVo partial 

K blackbody radiation. In that case, *^1“ »hown. Solid lines • 

the two curves corresponding to HIR case. Broken lines i LIR case, 

and L^R are the same. In the range 10^ S y^ i 5 x 10^ pair production losses 
play a role in decreasing y and therefore decreasing the fragmentation rate 
(u'le figure 1) so that the fragmentation time is somewhat longer and the curve 
is flatter. In the range y < 10^ pair production is again unimportant. 



u 


The combined effects of photodisintegration ^lnd pair production lead to 
attenuation of nuclei above a critical observed total energy as given for the 
six cases shown in table 1. Shown for comparison, is the critical energy for 
protons from photomeson losses (Stecker 1968) . 

TeUsle 1. Critical Attenuation Energy 



Univsrsal 

Origin 

(t.ixlO^^s) 

Supercluster Origin (t«10^®s) 

HIR 

4 

X lol® 

ev 

1 X 10^° eV 

LIR 

7 

X lOl® 

eV 

2 X 10^0 eV 

No IR 


X 10^^ 



(blac)dbody only) 

4 

eV 

no cutoff 

Protons 

(Photomeson losses) 

5.3 

X 10^^ 

eV 

no cutoff 



OU 
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For th« universal origin case, 
note that the presence of an 
Infrared photon field re- 

duces the critical energy by 
as much as an order of magnlt- 
vu'e below what would be expect- 
ed from the 2.7 K blackbody 
radiation alone. Furthermore, 
the numbers given In Table 1 
are derived for an initial com- 
position of pure Fe. For nuclei 
of lower Initial A , .tbe 
cutoff will be at a lower en- 
ergy. The predicted cutoffs 
for the universal origin case 
are very sharp. Since the data 
at present' are consistent with 
a smooth E“^ cosmic ray spec- 
trum up to 'V 1 or 2 X 10^^ eV, 
these data appear to rule out 
the universal origin hypothesis 
for bouh protons and nuclei but 
allow supercluster origin as 
well as galactic origin for 
protons and nuclei. 

Many more results have 
been obtained which will be 
discussed elsewhere (Puget, J.L., 
Stocker, F, W. and Brodekamp, 

J. J. , in preparation) . 
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Figure 2. Calculated history of a beam of 
Fe nuclei with Initial energy AMjjY for 
the HIR case for two values of y^. Also 
shown Is 1C where T_^ Is the life- 

time against one nucleon loss (see Figure 1) . 





Plgur* 3, Mean mass A(t) and spread In mass distribution 
(shaded area) for an initial Pe beam with three values of 
Yg. Solid lines I HZR case, Broken lines t LIR case. 
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